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Recent investigations in nonlinear sciences show that not only hyperbolic but also mixed dynamical systems 
may exhibit exponential relaxation in the chaotic regime. The relaxation rates, which lead the decay of proba- 
bility distributions and correlation functions, are related to the classical evolution resolvent (Perron-Frobenius 
operator) pole logarithm, the so called Pollicott-Ruelle resonances. In this Brief Report, the leading Pollicott- 
Ruelle resonances are calculated analytically for a general class of area-preserving maps. Besides the leading 
OO ■ resonances related to the diffusive modes of momentum dynamics (slow rate), we also calculate the leading 

f^**) ' faster rate, related to the angular correlations. The analytical results are compared to the existing results in the 

literature. 
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I. INTRODUCTION 



It is well known that for systems exhibiting chaotic dynamics, precise long-time predictions of individual trajectories are 
impossible. It is natural, therefore, to investigate the statistical properties of these systems. In this sense, the time evolution of 
^| ■ the probability densities of trajectories p„ , ruled by the Perron-Frobenius (PF) operator U as p n +i = Up n , have been extensively 
studied flU]. 

Due to Liouville's theorem, U can be represented by a unitary operator in a Hilbert space. Consequently, its resolvent 
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is singular on the unit circle in the complex z-plane, and the matrix elements of R(z) are discontinuous there. The sum in 
• (Q} is convergent for \z\ > 1 and has an analytical extension across the cut into the first Riemann sheet, which exhibits a set 
] of singularities known as Pollicott-Ruelle (PR) resonances JUlU]. To identify the PR resonances it is necessary to analytically 
f — ■ continue the resolvent across the continuous spectrum of U from the outside to the inside of the unitary circle. These resonances 
characterize the irreversible behavior of chaotic dynamics QUI]. In particular, the nontrivial (z ^ 1) maximal PR resonance 
. leads the exponential decay of distribution and correlation functions JalZl]- 

The PR resonances have attracted considerable attention not only in classical dynamics but also in quantum systems J^], 
C***'* , and some numerical and semi-analytical schemes where recently developed for calculate them. Blum and Agam proposed a 
variational method to locate the leading resonances JgD. Although their results describe the aparent formation of a leading 
\ quartet for two particular map cases, verified by respective numerical diagonalization of U, the leading resonance calculated 
diverges for a set of values of K in the standard map case when this approach breaks down. Florido et al. extended this 
variational approach in a class of numerical methods in which memory function and filter diagonalization techniques are utilized 
by means of interpolating exponentials ifioll . Usually, there are two standard ways to calculate the PR resonances: one is 
based on the numerical diagonalization of the operator U, for which the resonances are directly calculated from its eigenvalues 
llsil fllljl^l . the other, through the zeros of the classical Ruelle zeta-function, that is derived from the trace of the resolvent of 
U |Q1 0]. In the last case, there are analytical calculations of these resonances for some hyperbolic systems (for which this 
formalism is rigourous) such as the multibaker map [1], geodesic motion in billiards of constant negative curvature Ill3ll . and 
hard-disk scatterers [1]. On the other hand, many physically realistic systems are mixed, and analytical procedures to determine 
resonances for this cases are thus in demand. 

The motivation of the present Brief Report is to calculate analytically the leading PR resonances for slow (diffusive) and faster 
modes of dynamics for the general class of two-dimensional area-preserving maps: 

i n+1 = i n + Kf(e n ), 

9 n+ i = 9 n + ca(I n+1 ) mod27r, 
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defined on the cylinder — ir < 9 < tt, —oo < I < oo. Here f(9) is the impulse function, a{I) = a(I + 2irr) is the rotation 
number, c and r are real parameters, and K is the stochasticity parameter. This map is commonly called the radial twist map 
lfl4ll periodic in momentum variable /. The specific linear rotation number (LRN) case ca(I) = I for which f{9) = sm9 
represents the Chirikov-Taylor standard map 1131 . a paradigm of Hamiltonian chaos [ 14]. LRN maps are periodic because / can 
be replaced by / mod 2tt. On the other hand, nonperiodic rotation numbers can be considered in the limit r — ► oo 11611 . 



II. PROJECTION OPERATORS 



An usual way to determine the leading PR resonance is to evaluate the application U n for large values of the time n when 
only the highest resonance survives, as it occurs for the equilibrium statistical mechanics of lattice systems. Let us consider the 
analysis of the resolvent ([TJ for which U n can be expressed as § c dzR(z)z n = 2niU n J5[]. The spectrum of U is located on the 
unit circle C around the origin in complex z -plane or inside it. Thus, the contour of integration is a circle lying just outside the 
unit circle. In order to evaluate U n , a very effective method based on the projection operator techniques can be used 13, 1 1 811 . In 
this method, we consider two mutually orthogonal idempotent operators P and Q: 

1 = P + Q, P 2 = P, Q 2 = Q, PQ = QP = 0, (3) 

where 1 represents the identity operator. These operators decompose the resolvent in the following nontrivial form: 

' [P + QC{z)P] 1 — [P + PV(z)Q] + QV(z)Q, (4) 



z-U ^ w z-P£{z)P 

where the operators V{z), £{z), C(z) and V{z) are the discrete time version of the Brussels formalism [ 17], defined by 

QV(z)Q = Q i— Q, (5) 

z - QUQ 

P£{z)P = PUP + PUQV(z)QUP, (6) 

QC{z)P = QV(z)QUP, (7) 

PV{z)Q = PUQV(z)Q. (8) 

A recent proof of ©-([8]) can be found at ifTill . 

The matrix representation of the PF operator U for (f2]l in Fourier space (9, 1) t— * (m, q) is given by [7] 

(to, q\U\m\ l')=J2[ dq ' S ^ + q ^ Gl ^ mc )Jrn- m > {-Kq'), (9) 

m' I 

where the Fourier decompositions of the a(I) and f(9) functions are 

Ql(r,x) = J d9 cxp{-i[xa(r8) - 19}} , (10) 
Jm{x) = ^ J d9 G * V {-i[m9~xf{9)]}. (11) 



III. SLOW RELAXATION RATE 



The leading PR resonances related to diffusive modes of the momentum variable / for (O corresponds to the relaxation rate 
of PU n P ~ exp[nj(q)) for n ^> 1 and P = |0, q). The diffusion coefficient D is then calculated by D = — (l/2)d' 2 -f(q)\ q= o. 
Applying P on the two sides of the projection of the PF operator U n can be written as 

If z n 

PU n P = T1 f dz „oo (12) 
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where the memory functions &j(q) are given by [7]: 

MQ) = M-Kq), (13) 

*l(<z) = J2 J-m(-Kq)J m {-Kq)g (r,mc), (14) 
m 

*J>2(?) = X) S «^-tni(-- K '?)^n J (--«■«) 5Ai(r,mic) 
{m} {A}t 



X 
i=2 



fc=i 



(15) 



Hereafter, the following convention will be used: the set of wavenumbers to and {to} = {mi, . . . , rrij} can only take non-zero 
integer values, whereas the set of wavenumbers {A}^ can take all integer values, including zero, and the superscript denotes the 
constraint Yl&=i A* = 0. 

The integral dT2b can be solved by method of residues and its poles are evaluated by the well-known Newton-Raphson 
method: the zeros of an equation hjy(z) = are calculated iteratively by z„ +1 = z n — hjy(z n )/h' N (z n ), where h]\r(z) 
z — X)jLo z ^j(l) assumes the truncated form of the denominator of ( fT2b . First, we introduce the abbreviations M q = 
^2f =0 ^j(q) and N q = X)j=i 3 Notice that, taking into account the null drag condition J d6f(9) = 10], we have 

^o(q — > 0) = 1 + 0(q 2 ). In the general case, we have fyj>i(q — * 0) = 0(q 2 ). For q = 0, z* = 1 is the only root of h(z). 
This trivial pole is related to the equilibrium state found for to = m! = q = 0. For q — ► 0, the Newton-Raphson sequence of 
iterated roots will be given by zo = 1, z% = z<i = . . . = Zoo = M q + 0(q 4 ). For any choice of N > 1, it is easy to see that 
z 00 (N) is a root of the Hn(z), thus z* = limjv^oo z oa (N) is the leading pole of dl2l . Up to fourth order in q this pole can 
be considered simple because P£(z*)P = z* + 0(q 4 ). Performing the complex integration of (TT~2T > for n> lwe obtain the 
leading PR resonance 'y(q) H: 

oo 

1 (q)=lnJ2* ] (q)+0(q i ). (16) 

3=0 

The relaxation rate JToT ) is called slow because 7(g) = 0(q 2 ) for small wave number q. 

IV. FASTER RELAXATION RATE 

Likewise the leading resonance corresponding to the diffusive modes of the momentum variable / leads the exponential relax- 
ation of distribution functions, leading angular resonances have a important role in the exponential decay of angular correlation 
functions 

C uv (n) = (u\U n \v) ~ e~ n \ (17) 

where u and v are two any observables at the same instant of time and n is sufficiently large. Let us consider the analysis of the 
transition elements Q\U n Q = (m, 0\U n \m', q'). The expansion of QiR(z)Q can be written as 

Qi — rr Q = V^ (4+1) ^, (18) 

Z — U *7~t 

2 — 1 

where fa = Q\XJ l Q. The analysis becomes simpler for the LRN case, where Q\(r = 1, x) = S\ jX , and for which we have the 
following first three fa coefficients: 

fa = ^Jtn-m'i-mK), (19) 
m' 

fa = ^2j m -x(-mK)^2jx- m '[-(m + ^)K], (20) 

A m' 

fa = J2 m (-mK)Y,J-( m + m ')(mK) + T m {K) + 0{J 3 ), (21) 
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where 

T m {K) = Jl(-mK) + J (-mK)J m (-mK) 

+ Jm-mi-mK) {J m+2m > [- (m + m')K] + J m , [-(m + m')K}} . (22) 

m' 

In the calculation of ([T9ll-(l2"Tli. as well as in the calculations that follow, it is crucial to consider the following addition rule 

y.Jm-m'jx) = J^Jm-xjx) - J m (x) = 1- J m (x). (23) 

rn' X 

Notice that, including I = 0, we have J^i exp(— iif) — 2ttJ2i ^(t ~ 2ttI). Hence, the identity d23l holds due to ^ A J\{x) = 
exp[ia:/(0)] = 1 for /(0) = 0. Such a result was only known for the particular case of Bessel functions of first kind by means 
of its generating function 1 19]. 

For sufficiently high values of K we expect that the coefficients fa become negligible as i increases. Thus, in a first approxi- 
mation, we can truncate the right hand side of ( fT8l at i = 3 and rewrite it in the following rational form 



M 2 + 022 



(24) 



whose coefficients ipi are given in terms of fa as 



; 1 ; 02 , 02 91 

03 03 03 03 



The right hand side of d24b is, in a first approximation, the analytical extension of the series representation of Q\R(z)Q, valid 
in the chaotic regime. The non null poles of the projected resolvent d24"l i form the leading resonances of the PF operator Q\UQ. 
First, we have fa = 1 — J m (—mK) due to d23l . thus fa ^ unless K = — 1 for the particular case of the sawtooth map 
f(6) = 9. Considering fa = 1 + O(J) as the dominant term, fa (for A = m' = — rri) and fa must be the O(J) perturbative 
terms of the (^-expansion d24"l i. Neglecting only (D(J 3 ) terms on the fa coefficients, the poles of the rational form (l24i i will be 
given by 

z± = ±M-ltl + 0(j3/>). (26) 
V 01 2 fa 

The ratio fa/ fa can be considered only as J^2m(—mK), and its 0(J 2 ) terms can be neglected. On the other hand, the ratio 
fa/ fa ~ ■J2m(—'mK) must be considered up to <D(J 2 ) terms, given C( l / 1 / 2 ) and 0{J) corrections, respectively. Thus, the 
leading angular resonance, represented in the exponential form as \z\ = cxp(— 7), will be 



In max 



1 



y/j2m(-mK) + T m (K)± -J 2m (-mK) 



(27) 



with T m (K) given by (1221 . Note that, for odd impulse function f(6), the leading resonance d27l i is invariant under the change 

m — > —m. 

For very large values of K, the leading resonance ( |27] | tends to the following value: 

7oc = - lnmax\y/ J 2m (-mK)\, (28) 

rn 

obtained in a different way by Khodas et al. for the standard map particular case |@] . It is important to check the limits of validity 
of each approximation and its respective adequacy to the numerical values existing in the literature. In Fig. Q]we compare, for 
the standard map, the resonance d27l i with its asymptotic value (l28l . For sufficiently large times, Khodas et al. J6J calculate 
numerically the correlation (fTTT i for u and v proportional to exp(im^) and for some combinations of modes {m, m'}, where 
C U v = C m ,m' in this choice. Once the resonance ( |27] i is dominated by initial mode rn = 1, we select C\ \ and C1.2 as the 
best simulated correlations. However, these numerical values have only a qualitative character for sake of comparison, since 
resonance d27] > leads the decay of correlations only for very large times, when the numerical signal is too weak [6]. Moreover, 
the precise composition of the observable v as a possible superposition of modes {m^} for which C\ jV decay through d27l > is not 
known. On the other hand, we also include the two values of leading resonances calculated numerically by diagonalization of U 
for K — 10 and K — 13 by Blum and Agam [9], in addition to the leading intermediary value calculated by Florido et al. for 
K = 10 (assumed here between z — 0.672 and z = 0.715) floll . By comparing all these results, the theoretical result ( l27b gives 
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FIG. 1: Theoretical leading resonance 7 (solid line) calculated for the standard map compared with its asymptotic value 700 (dotted line) 
and several numerical calculations. Here, (■) and (A) represent the resonances calculated from Ci,i and Ci,2 correlations, respectively, by 
Khodas et al. @], (*) represents the resonances calculate by Blum and Agam [9], and (if) is the intermediary value calculated by Florido et 
al. (U. 



a better qualitative fit with the numerical values even for high values of K and, besides, it reveals a more intrincated structure of 
peaks for the maximal resonance. 

The sawtooth map f(9) = 9 is the only particular LRN case for which the perturbative scheme presented above breaks down 
02OI 1. This occurs in such case due to J m (x) = sin J^- x)1 = 1 for m = x and interger values of K. For example, besides 
4>i = 1 for all integer K ^ — 1, <f>2 can be rewritten as 

fo = l-J2Jrn-\(-mK)J x [~(m + \)K}. (29) 

A 

For integer K, the sum in the right hand side of (l29t vanishes unless —K= m ^ A = m m X , which gives |A/m| = g» or 

|A/m| = g^ 1 , where = (\/5 — l)/2 is the golden mean. Hence, we also have (f>2 — 1 for all integer K. This suggests that 
<f>i = 1 for almost all integer K. If this hypothesis is true, we then have as analytical continuation of the resolvent for \z\ > 1 

^lhj Q = V^ = Khry (30) 

3=0 V ' 

according to ( fl~8T >. Thus, z = 1 corresponds to the invariant density and all the others resonances are infinitely degenerated at 
2 = 0. This particular result was demonstrated by Sano for all positive integers K by using the Fredholm determinat of U IU2I1 . 



V. CONCLUDING REMARKS 



In conclusion, we have presented a method to determine analytically leading Pollicott-Ruelle resonances which is applicable 
to a general class of area-preserving maps, including mixed systems. Such resonances are obtained through the resolvent of 
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the PF operator by using projection operator techniques. In particular, we calculate the leading resonance related to the slow 
modes of relaxation, which corresponds to the diffusive process, as well as the leading resonance related to the faster modes of 
relaxation. In this last case, our perturbative analysis was performed only for systems with linear rotation numbers, although it 
can be similarly applied for nonlinear ones. 

The analytical results obtained here have been compared with theoretical and numerical calculations existing in the literature. 
The resonance d27b was calculated in a systematic way in which correction terms of order O(J) produce a more intrincated 
structure of peaks for the standard map case even for high values of K, as can be seen in Fig. Q] Despite the absence of estimates 
of errors in the numerical results, the agreement with the theoretical result d27l i is reasonable. We have also investigated particular 
characteristics of the sawtooth map that are incompatible with the perturbative approach developed in the Sec. [IV] Our analysis 
points towards the accordance between our hypotesis and the results presented in 11211 . 
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